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Scaling behavior of the conserved transfer threshold process 

S. LiibeckQ 

Weizmann Institute of Science, Department of Physics of Complex Systems, 76100 Rehovot, Israel, 
Theoretische Tieftemperaturphysik, Gerhard-M creator- Universitdt, 47048 Duisburg, Germany 

(Dated: June 20, 2002) 

We analyze numerically the critical behavior of an absorbing phase transition in the conserved 
transfer threshold process. We determined the steady state scaling behavior of the order parameter 
as a function of both, the control parameter and an external field, conjugated to the order parameter. 
The external field is realized as a spontaneous creation of active particles which drives the system 
away from criticality. The obtained results yields that the conserved transfers threshold process 
belongs to the universality class of absorbing phase transitions in a conserved field. 
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I. INTRODUCTION 

The scaling behavior of directed percolation is recog- 
nized as the paradigm of the critical behavior of sev- 
eral non-equilibrium systems which exhibits a continu- 
ous phase transition from an active state to an absorbing 
non-active state (see for instance The widespread 

occurrence of such systems in physics, biology, as well as 
catalytic chemical reactions is reflected by the well known 
universality hypothesis of Janssen and Grassberger that 
models which exhibit a continuous phase transition to a 
single absorbing state generally belong to the universality 
class of directed percolation |2[ || . Introducing additional 
symmetries the critical behavior differs from directed per- 
colation. In particular particle conservation leads to the 
different universality class of absorbing phase transitions 
with a conserved field as pointed out in Q . In that work 
the authors introduced two models, the conserved lattice 
gas (CLG) and the conserved threshold transfer process 
(CTTP). The latter one is a conserved modification of the 
threshold transfer process introduced in Q . Both models 
display a continuous phase transition from an active to an 
inactive phase and are believed to belong to the same uni- 
versality class The steady-state scaling behavior of 
the CLG model was investigated recently. The order pa- 
rameter and its fluctuations were numerically examined 
in The scaling behavior in an external field conju- 
gated to the order parameter was considered in [Q. Fur- 
thermore a modified CLG model was introduced which 
allows to determine analytically the steady-state mean- 
field scaling behavior of the universality class j|, ||] . 

On the other hand the scaling behavior of the CTTP 
was investigated in low dimensional {D ~ 1,2) systems 
only |l0| and no external field was applied. Therefore 
we consider in this work the CTTP with and without an 
external field in various dimensions {D = 1,2,3,4,5,6) 
and determine a set of critical exponents. All obtained 
results coincides with those of the CLG model, strongly 
supporting the universality hypothesis of B . 



II. MODEL AND SCALING BEHAVIOR 

In this work we consider the CTTP on simple cubic 
lattices of linear size L in various dimensions with N 
particles. The lattice sites may be empty, occupied by 
one particle, or occupied by two particles. Empty and 
single occupied sites are considered as non-active whereas 
double occupied lattice sites are considered as active. In 
the latter case one tries to transfer both particles of each 
active site to randomly chosen empty or single occupied 
nearest neighbor sites. If no active sites exist the system 
is trapped forever in a certain configuration, a so-called 
absorbing state. 

In the following we denote the densities of active sites 
with Pa and the density of particles on the lattice as 
p = N/L-^, which is considered as the control param- 
eter of the absorbing phase transition. The density of 
active sites is the order parameter of the absorbing 
phase transition, i.e., it vanishes at the critical density 
Pc according to 



Pa ^ 6p^, 



(1) 
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with the reduced control parameter Sp — p/ Pc — 1 and 
the order parameter exponent f3. 

Similar to equilibrium phase transitions it is possible 
in the case of absorbing phase transitions to apply an ex- 
ternal field h which is conjugated to the order parameter 
(see for instance Q). As usually for continuous phase 
transitions the conjugated field has to destroy the dis- 
ordered phase and the associated linear response func- 
tion dpa/dh has to diverge at the critical point {Sp = 0, 
/i = 0). In the case of an absorbing phase transition the 
external field acts as a spontaneous creation of active par- 
ticles, i.e., the external field destroys the absorbing state 
and thus the phase transition itself. But considering ab- 
sorbing phase transitions with particle conversation one 
has to take care that the external field does not change 
the particle number. A possible realization of the ex- 
ternal field was developed in [Q where the external field 
triggers movements of inactive particles which may be 
activated in this way. The external field h is another rel- 
evant scaling field and for sufficiently small values of h 



2 



the order parameter scales as 

p^{5p,h) ^ Xf{SpX~^^^,hX--/^) (2) 

with the critical field exponent a and the scahng func- 
tion r. Choosing SpX^^^^ = 1 one recovers Eq. (|l|) 
whereas h X~"/^^ = 1 leads at the critical density to 



Pa 



(3) 



In our simulations we start with randomly distributed 
particles. All active sites are listed and this list is up- 
dated in a randomly chosen sequence. In the case that 
an external field is applied the active particle creation 
is performed after each update step in order to mimic 
the external field. After a certain relaxation time the 
system reaches a steady state where the density of active 
sites at update step t fiuctuates around the average value 
{pa_{5p^h,t)) which is interpreted as the order parameter 
Pa{5p, h) (see for instance Figs. 1 of 0]). 

Additionally to the order parameter we consider its 
fiuctuations 

l^p,{5p,h) = [{p,{6p,h,tf) - {p^{5p,h,t)f]. 

. . (4) 

Approaching the transition point the fluctuations diverge 
for zero-field according to 



l\p^{5p,h = {)) - 5p- 



(5) 



The fluctuation exponent 7' fulfills the scaling rela- 
tion ill 



i - v^D - 2/3, 



(6) 



where the exponent Vj^ describes how the spatial correla- 
tion length diverges at the transition point. In the critical 
regime we assume that the fiuctuations obey the scaling 
ansatz 



Ap^{Sp,h) = A'' d{SpX,hX''). 



(7) 



Setting 5pX—l one recovers Eq. (||) for h — 0. 

Analogous to equilibrium phase transitions the suscep- 
tibility is defined as the derivative of the order parameter 
with respect to the conjugated field 



d 

x{Sp,h) = —p^{5p,h) 



= X^--ll' diSpX-^'P.hX-'^'^). 



(8) 



Setting 5pX ^1^^ — 1 one gets that the susceptibility di- 
verges for zero-field as required according to 



dh 



Sp- 



Furthermore, one yields the scaling relation 
7 = cr - /3 



(9) 



(10) 



which corresponds to the well known Widom equation of 
equilibrium phase transitions. Using this scaling relation 
one can calculate the value of the susceptibility expo- 
nent 7 from the obtained values of (3 and a. Notice that 
in contrast to the scaling behavior of equilibrium phase 
transitions the non-equilibrium absorbing phase transi- 
tion is characterized by 7 7^ 7'. 



III. BELOW THE UPPER CRITICAL 
DIMENSION 

At the beginning of our analysis we consider the scaling 
behavior of the order parameter for D = 1, 2, 3. System 
sizes up to i = 131072 for D = 1, L = 2048 for D = 2, 
and L = 256 for {D = 3) are considered. In each cases 
we start the simulation with randomly distributed parti- 
cles. After a certain transient regime the system reaches 
a steady state where the density of active particles fluctu- 
ates around an average value which is interpreted as the 
order parameter. In the steady state up to 2 10^ update 
steps for D = 1 and 2 10^ for D = 2, 3 are performed to 
measure the average density of active sites. For zero-field 
this procedure is repeated for at least 10 different initial 
configurations in order to get an accurate estimation of 
the order parameter close to the critical point {p — pc, 
h = 0). 

In Fig.|l| we present the data of the one-dimensional 
order parameter at zero-field. Approaching the transition 
point the corresponding correlation length increases and 
the system tends to the absorbing state if the correlation 
length is of the order of the system size. Instead of a 
finite-size scaling analysis (see for instance |l^, |2)) 
we take care of these finite-size effects in the way that 
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FIG. 1; The order parameter pa as a function of the par- 
ticle density for zero-field (symbols) and for h — 0.0001 and 
h — 0.00005 (lines). The inset displays the order parameter 
fluctuations Apa for zero field (symbols) and for h — 0.0002, 
h = 0.0001 and h = 0.00005 (lines). 
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FIG. 2: The order parameter pa as a function of the reduced 
particle density 8p at zero-field for various dimensions D. The 
dashed line corresponds to a power-law behavior according to 
Eq. (|^) for D ^ D^. For D — Q the data are shifted horizon- 
tally by a factor 1.5 in order to avoid an overlap. In the case 
of the four-dimensional model the dashed line corresponds to 
Eq. with B = 0.15. 



we increase the system size before these finite-size effects 
occur and use only data from simulations that have not 
reached the absorbing state. 

Decreasing the particle density the order parameter 
decreases and vanishes at the transition point. To de- 
termine the critical indices one varies the critical den- 
sity pc until one obtains asymptotically a straight line in 
a log-log plot. The exponent is then obtained by a re- 
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FIG. 3: The order parameter fluctuations Apa as a function 
of the reduced particle density dp at zero-field for various 
dimensions D. The dashed line corresponds to a power-law 
divergence [Eq. (^]. For D > Dc the fluctuations are maximal 
at the transition point but finite. 



FIG. 4: The scaling plot of the order parameter fluctua- 
tions Apa for various dimensions. For D > 1 the curves 
are shifted vertically in order to avoid overlaps. In the case 
of the four-dimensional model pa ft"^''^ | In is plotted 
vs. 5ph~^^^\lnh\''-''^'^ with S = 0.28 and 6 - s/2 = -0.12 
(see text). 



gression analysis. The values of the order parameter as 
a function of the reduced particle density Sp are plotted 
m Fig. I In all cases the asymptotic behavior {Sp 0) 
of the order parameter obeys Eq. (|^) . For Z? = 1 we get 
Pc = 0.96929 ± 0.00003 and (3 = 0.382 ± 0.019. The lat- 
ter value is smaller than the value (3 — 0.412 estimated 
from significantly smaller system sizes {L < 5000) [ p^ . 
Furthermore our value differs from (3 — 0.42 ± 0.02 ob- 
tained from simulations of the one-dimensional fix-energy 
Manna sandpile model |l^ that is expected to belong to 
the same universality class. 

In the two dimensional case we obtain pc — 0.69392 ± 
0.00001 and /3 = 0.639 ± 0.009. Again the order pa- 
rameter exponent differs slightly from the previously re- 
ported result P — 0.656 obtained from simulations of 
small lattice sizes {L < 160) |jl^. But our value agrees 
with the estimate of the corresponding two-dimensional 
Manna sandpile model P = 0.64 ± 0.01 ||l||. 

The estimates of the three dimensional model are pc = 
0.60489 ± 0.00002 and /3 = 0.840 ± 0.012. All obtained 
critical exponents are listed in Table ||. 

In Fig. ^ we present the order parameter fiuctuations 
as a function of the control parameter at zero field. We 
observe for D < Dc a power-law behavior according to 
Eq. (||). Using a regression analysis we get the estimates 
Y = 0.662 ± 0.071 for £) = 1, 7' = 0.381 ± 0.013 for 
£1 = 2, and 7' = 0.187 ± 0.030 for D = 3. 

In the following we analyze the order parameter as a 
function of the control parameter Sp for different fields 
from h — 10~^ up to 2 10~^. The applied field results 
in a smoothing of the zero-field curve, i.e., the order pa- 
rameter increases smoothly with the control parameter 
for h > (see Fig.|l|). According to the scaling ansatz of 
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FIG. 5: The scaling plots of the order parameter fluctua- 
tions AjOa for various dimensions. For D > 1 the curves are 
shifted vertically in order to avoid overlaps. The fluctuations 
diverges at the critical point for D < 4 whereas a jump of 
the fluctuations is observed in higher dimensions at zero-fleld. 
In the in case of the four- dimensional model Apa is plotted 
vs. Jp/i"^/^ lln/ip" with T] = 0.1. 



the order parameter [Eq. 
get the scahng form 



we choose h A '^^^ — 1 and 



(11) 



Thus one varies the exponent a until the curves for dif- 
ferent values of the driving field have to collapse onto 
the scaling function f if one plots pa h~^^'^ as a func- 
tion of 5ph~^/'^. Convincing results are obtained for 
a = 1.770 ±0.058 {D = 1), ct = 2.229 ±0.032 {D = 2), as 
well as (7 = 2.069 ±0.043 [D — 3) and the corresponding 
scaling plots are shown in Fig.|^. 

Next we consider the scaling behavior of the order 
parameter fluctuations Apa- The fluctuation data for 
D — \ are shown for different values of the external 
field in the inset of Fig.|^. For finite fields the fiuctu- 
ations display a peak. Approaching the transition point 
{h —>■ 0) this peak becomes a divergence signalling the 
critical point. In order to analyze the scaling behavior of 
the fluctuations we use the scaling ansatz Eq. and set 



Api,{Sp, h) 



d{5p h 



-l/a 



I) 



(12) 



Using the above determined values of pc, cr and 7' we 
get good data collapses conflrming the accuracy of our 
analysis, (see Fig.^) 

Furthermore we determine the susceptibility expo- 
nent 7. Using the scaling relation Eq. ([l^) one gets the es- 
timates of the susceptibility exponents 7 — 1.388 ±0.063 
(£) = 1), 7 = 1.590±0.033 (£» = 2), and 7 = 1.229±0.045 
(7? = 3). 



IV. AT THE UPPER CRITICAL DIMENSION 

In the case of the four dimensional model we considered 
system sizes from L = 8 up to L = 64. At least 10^ up- 
date steps were used to reach the steady state close to the 
transition point and 2 10^ update steps were performed 
to determine the order parameter and its fluctuations. 

At the upper critical dimension Dc = 4 the scaling 
behavior of the CTTP is affected by logarithmic correc- 
tions similar to the CLG model H, 0. As argued in |^ 
the order parameter obeys in leading order the scaling 
ansatz 



Pi,{Sp,h) = 



(13) 



A I lnA|' f{SpX-^/f^\ lnA|^ A-'^/'^l lnA|^), 



where the exponents /3 and a are given by the correspond- 
ing mean-fleld values /3 — 1 and cr = 2 respectively. 
Thus, for zero field the asymptotic scaling behavior of 
the order parameter obeys 



p^{Sp,h^O) ~ 6p\\n6p\^ 



(14) 



with B = b+l. In our analysis we plot pa/ Sp as a function 
of I In Zip and vary the exponent B as well as the critical 
density pc until one gets asymptotically a straight line. 
The best resuh is obtained for B = 0.15, pc = 0.56705 ± 
0.00003 and the corresponding plot is shown in Fig.||. 
This value of B differs from the corresponding value of 
the CLG model B = 0.24 Q. 

Similar to the lower dimensions we consider the scal- 
ing behavior of the order parameter as a function of the 
control parameter for different external fields. Choosing 
h\-''/l^\\nX\'' = 1 the scaling ansatz [Eq. (|l|)] yields in 
leading order 

(15) 



0.60 



0.55 



10 

CO 

Q. 



0.50 



0.45 



Pa('5p, h) 


= |ln/i|^ f(a;,l). 


. L=8 




■ L=16 
♦ L=32 


y 

X 

y 


A L=64 


y 

* 

y 






y 

y 

y 

y 

y 


y 

y 

y 

D=4 



1.00 



1.10 



|ln5p|' 



1.30 



FIG. 6: The density of active sites at the upper critical di- 
mension Dc = 4. The data are rescaled according to Eq. (p^. 
The assumed asymptotic scaling behavior (dashed line) is ob- 
tained for B = 0.15. 
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where the scahng argument x is given in leading order by 

X = 5ph-'^/'^\\iih\^-'/'^ (16) 

with S — s/2 + I. Varying the logarithmic correction 
exponents one gets for S = 0.28 and b — s/2 = —0.12 a 
convincing data-collapse, which is shown in Fig.Q. Using 
the values E = I + s/2 = 0.28 and b - s/2 = -0.12 
we get the estimation B = b + I — 0.16 which agrees 
with B — 0.15 obtained from numerical simulations in 
zero-field. On the other hand this values differs from the 
corresponding estimations of the CLG model B = 0.24, 
E = 0.45, and 6 - s/2 = -0.17 0. 

Furthermore we consider how the logarithmic correc- 
tions affect the scaling behavior of the fluctuations at the 
upper critical dimension. As pointed out in j^] the order 
parameter fluctuations are expected to obey the scaling 
ansatz 



10" 



Apa{5p,h) = d{5ph 



-1/2 



(17) 



A good data collapse is observed for ?/ — 0.10 (see Fig.|^) 
which differs again from the corresponding value of the 
CLG model tj = 0.39 §. 



V. ABOVE THE UPPER CRITICAL 
DIMENSION 

A modified version of the CTTP with random neigh- 
bor hopping was recently introduced in There, unre- 
stricted particle hopping breaks long range correlations 
and the scaling behavior is characterized by the mean- 
field values pc — 1/2, (3 — 1, and (7 = 2 which are calcu- 
lated analytically. 

In our simulations of the five dimensional model we 
considered system sizes from L 8 up to L 32 whereas 
system sizes from L = 4uptoL = 16 are used for D = 6. 
At least 2 10^ update steps were used to reach the steady 
state and 2 10^ update steps were performed to determine 
the order parameter and its fluctuations. The values of 
the order parameter are plotted in Fig. ^ and the obtained 
critical densities are pc = 0.54864±0. 00005 for D = 5 and 
Pc = 0.53816 ± 0.00007 for D = 6, respectively. In both 
dimensions the asymptotic scaling behavior of the order 
parameter is in agreement with the mean-field behavior 
(3 = 1. 

The fluctuations of the order parameter Apa are plot- 
ted in Fig.^. Analogous to the CLG model the fluctua- 
tions are characterized by a jump at the transition point 
corresponding to 7' = Q. 

Above the critical dimension, i.e. -D > 5, the scaling 
behavior of the CTTP is expected to obey again the scal- 
ing ansatzes Eqs. (^,0) where the exponents are given by 
the mean-field values independently of the particular di- 
mension. The obtained data collapse of the order param- 
eter curves are presented in Fig. ^ and confirm the above 
scenario. 

Furthermore we consider the fluctuations above the up- 
per critical dimension. According to the mean-field value 
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FIG. 7: The critical density pc{D) as a function of the di- 
mension D. The critical density of the mean-field solution is 
denoted by Pc(cxd) = 1/2. The dashed line corresponds to a 
power-law behavior [Eq. with an exponent 1.48. 



7' = we plot Apa as a function of 6ph^^/^ and the 
obtained data collapses are shown in Fig. ||. 

Finally we address the question how the critical densi- 
ties depends on the dimension. As can be seen from Ta- 
ble|| the critical density tends with increasing dimension 
to the mean-field value pc = 1/2 ||] that corresponds to 
an infinite dimension. Our analysis reveals that the crit- 
ical densities approaches that mean-field value according 
to 



D- 



(18) 



with T = 1.48 ± 0.05 (see Fig.0). This behavior is dif- 
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FIG. 8: The critical exponents 0, 7', and cr of the CTTP and 
CLG model (obtained from ^,0) for various dimensions. The 
dashed lines are just to guide the eyes. 
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ferent from that of CLG models on simple cubic lattices 
which is characterized by an exponent r = 1 (§1 . 



VI. CONCLUSIONS 

We investigated the steady state scaling behavior of 
the CTTP model in various dimensions. The order pa- 
rameter exponent, the fluctuation exponent and the ex- 
ternal field exponent are determined and the correspond- 
ing values are Hsted in Table |. For D = 1 and Z? 2 
our results of the order parameter exponents differ from 
previous simulations obtained from significantly smaller 
system sizes. Our values of the critical exponents (3, 7', 
and a agree within the error-bars with the corresponding 
exponents of the CLG model (see Fig.^, strongly sup- 
porting the conjecture Q that both models belong to the 
same universality class. 

The picture is not so clear at the upper critical di- 
mension Dc = 4. Although the exponents are identical 
the logarithmic correction exponents of the CTTP and 
CLG model are different. This result is rather surprising 
since the logarithmic corrections exponents are a char- 
acteristic feature of the whole universality class (see for 
instance 0). We think that more than statistical un- 
certainties this result is caused by systematic uncertain- 
ties of our analysis. In all cases we focused our atten- 
tion to the leading order of the scaling behavior. Taking 



corrections to the leading order into account may result 
in comparable values of the logarithmic correction expo- 
nents. Further investigations are needed to clarify this 
point. 
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TABLE I: The critical density pc and the critical exponents 
/3, a, 7' and 7 of the CTTP model for various dimensions D. 
The values of the susceptibility exponent 7 are calculated via 
Eq. ( |lo| ) . The symbol * denotes logarithmic corrections to the 
power-law behavior. 
D Pc P (T 7^ 7 



1 0.96929 0.382 ± 0.019 1.770 ± 0.058 0.662 ± 0.071 1.388 ± 0.063 

2 0-69392 0.639 ± 0.009 2.229 ± 0.032 0.381 ± 0.013 1.590 ± 0.033 

3 0-60489 0.840 ± 0.012 2.069 ± 0.043 0.187 ± 0.030 1.229 ±0.045 

4 0-56705 1* 2* 0* 1* 

5 0-54864 1 2 1 

6 0-53816 1 2 1 
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